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Abstract 

We present the first observation of quantum interference in the process (p^KsKi 
-^7r"'"7r~7r+7r~, using the KLOE detector at the Frascati e~^e~ collider DA$NE. 
From about 5 x 10^ neutral kaon pairs both decaying to 7r"'"7r~ pairs we obtain the 
distribution of At, the difference between the two kaon decay times, which allows 
testing the validity of quantum mechanics and CPT invariance: no violation of 
either is observed. New or improved limits on coherence loss and CPT violation are 
presented. 



1 Introduction 

A 0-factory provides unique opportunities for testing quantum mechanics (QM) 
and CPT symmetry. In the decay (p —>■ K^K^, the neutral kaon pair is produced 
m a J^^ = 1— state: 
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where p is the kaon momentum in the meson rest frame, and A^= (l + |ep)/(l — 
e^). Since 1 — iV <C 1 we will set iV = 1 in the following without any loss 
of generality. The decay intensity for the process 0— *-(2 neutral kaons)— ^tt+tt", 
n~^7i~ is then given by [1]: 
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where ti are the proper times of the two kaon decays, Ts and Tl are the de- 
cay widths of Ks and K^, Am = rriL — ms is their mass difference and ?7+_ = 
{7r~^7r~\Ki) / {7r~^ii~\Ks) = |?7+_|e*'^+". The two kaons cannot decay into the same 
final state at the same time, even though the two decays are space-like sepa- 
rated events. Correlations of this type in QM were first pointed out by Einstein, 
Podolsky, and Rosen (EPR) [2]. 

While it is not obvious what a deviation from QM might be, the assumption 
that coherence is lost during the states time evolution does violate QM. One 
can therefore introduce a decoherence parameter ( [3], simply multiplying the 
interference term in Eq. (2) by a factor of (1 — C). The meaning and value of 
( depends on the basis in which the initial state (1) is written [4]. Eq. (2) is 
modified as follows: 
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in the K^-K"^ basis. 



Another phenomenological model [5] introduces decoherence via a dissipative 
term in the Liouville-von Neumann equation for the density matrix of the state 
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and predicts decoherence to become stronger with increasing distance between 
the two kaons. This model introduces a parameter A, related to the decoherence 
parameter in the Ks-Kl basis by Csl — A/F^. 

In a hypothetical quantum gravity, space-time fluctuations at the Planck scale 
(~10~^^ cm), might induce a pure state to become mixed [6]. This results in QM 
and CPT violation, changing therefore the decay time distribution of the K^- 
pair from (p decays [7]. Three CPT- and QM- violating real parameters, with 
dimensions of mass, a, (3, and 7, are introduced in [7]. and 7 are guessed to 
be of 0{mj^/Mp) ~ 2 x 10"20 GeV [7,8], where Mp = = 1.22 x lO^*' GeV 

is the Planck mass. The conditions = 7 and f3 = ensure complete positivity 
in this framework [9-11]. The decay intensity is (see Eq. (7.5) in Ref. [9] setting 
a = 7 and /? = 0): 

1 , ^ , 4 , ,2 



[ (1 + j {^-''''-'''' + e-^^*-^^'^) (5) 
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It has been pointed out [12, 13] that in this context the initial state (1) may 
acquire a small C-even component: 

\^) = ^ (\K\ p)\K^-p) - \K', p)\K', -p) 



V2 

+ uj{\K', p)\K',-p) + \K', p)\K', -p))) (6) 
where u = |ci;|e*^ is a complex parameter describing CPT violation, whose order 



of magnitude is expected to be at most \uj\ ~ ^{mj^/Mp)/ AT ~ 10 ^, with 
AF = F5 — F^. The decay intensity is (see Eq.(3.3) in Ref. [13] setting a, 7 = 0): 
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The decoherence parameters (sl and Coo; have been found in the past to be com- 
patible with zero, with uncertainties of 0.16 and 0.7, respectively, using CPLEAR 
data [4,5,14]. CPLEAR has also analyzed single neutral- kaon decays to measure 
the a, P, and 7 parameters [15]. The values obtained for all three parameters are 
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compatible with zero, with uncertainties of 2.8 x 10"^^ GeV, 2.3 x 10"^^ GeV, 
and 2.5 x 10"^^ GeV, respectively. The parameter uj has never been measured. 

In the following, the improved KLOE measurements of the Csl, Coo; -^j 7; ^"^^ ^ 
parameters are presented. The analysis is based on data collected at DA$NE in 
2001-2002, corresponding to an integrated luminosity of L ^ 380 pb~^. DA$NE, 
the Frascati factory, is an e^e~ collider operated at a center of mass energy 
W = M(</))~1020 MeV. Electrons and positrons collide in the horizontal plane 
at an angle of tt— 25 mrad. 0-mesons are produced with approximately 12 MeV/c 
momentum toward the rings center, along the x-axis. The 2;-axis is taken as the 
bisectrix of the two beams, the y-axis being vertical. 

2 The KLOE detector 

The KLOE detector consists of a large, cylindrical drift chamber (DC), sur- 
rounded by a lead/scintillating- fiber electromagnetic calorimeter (EMC). A super- 
conducting coil around the calorimeter provides a 0.52 T field. The drift cham- 
ber [16] is 4 m in diameter and 3.3 m in length. The momentum resolution is 
Cpx/P-L ~ 0.4%. Two-track vertices are reconstructed with a spatial resolution 
of ~ 3 mm. The calorimeter [17] is divided into a barrel and two endcaps. It 
covers 98% of the solid angle. Cells close in time and space are grouped into 
calorimeter clusters. The energy and time resolutions for photons of energy E 
are ds/E = 5.7%/y^7GeV) and at = 57 ps/ ^JE~iGeV) ® 100 ps, respectively. 
The KLOE trigger [18] uses calorimeter and chamber information. For this anal- 
ysis, only the calorimeter signals are used. Two energy deposits above threshold 
{E > 50 MeV for the barrel and > 150 MeV for the endcaps) are required. 

Kaon regeneration in the beam pipe is a non negligible disturbance. The beam 
pipe is spherical around the interaction point, with a radius of 10 cm. The walls 
of the beam pipe, 500 fim thick, are made of a 62%-beryllium/38%-aluminum 
alloy (AlBeMet®162). A beryllium cylindrical tube of 4.4 cm radius and 50/im 
thick, coaxial with the beam, provides electrical continuity. 

We only use runs satisfying basic quality criteria. For each run we determine 
the average collision conditions: p=Pe++Pe =P4>, the center of mass energy W, 
the beam bunch dimensions, the collision point rc and angle. This is done using 
Bhabha scattering events. We then require \py,z\ < 3 MeV, \W — 1020 1 < 5 MeV. 
The collision point must satisfy | | < 3 cm, | yc I < 3 cm, and \zc\ < 5 cm. The 
rms spread of the luminous region must satisfy (Xx < 3 cm and az < 3 cm. A small 
number of runs were rejected because of improper trigger operation. Each run 
used in the analysis is simulated with the KLOE Monte Carlo (MC) simulation 
program, GEANFI [19], using values of relevant machine parameters such as W 
and determined as mentioned above. Machine background obtained from data 
is superimposed on MC events on a run-by-run basis. For KsKl^'^~^t^~'^~^'^~ 
events, the number of simulated events is 10 times that expected on the basis of 
the integrated luminosity. For all other processes, the effective statistics of the 
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simulated sample and of the data sample are approximately equal. The effects 
of initial- and final-state radiation are included in the simulation. Final-state 
radiation in Ks and Kl decays is treated as discussed in Ref. [20]. 

3 Analysis 

3.1 Event selection 

Events are selected on the basis of two identified neutral kaons from </> decay, in 
turn decaying into tt+tt" pairs. Since we cannot tell whether tt+tt" decays near 
the interaction point (IP) are from a Ks or a Kl, we try to ensure that the same 
criteria are used for all decays. In the following we call Ki {K2) the decay closest 
to (farthest from) the (j) production point. 

We first require a Ki vertex with two tracks of opposite curvature within a fiducial 
volume with r < 10 cm and \z \ < 20 cm centered at the nominal collision point 
rc, determined as discussed above. We also require that the two tracks satisfy 
I m^+TT- — niK I < 5 MeV and | p^+^r- — Pk \ < MeV/c, where px is calculated 
from the kinematics of (p^KsK^. m^+^- and respectively the invariant 

mass and the momentum of the tt^ti~ pair. 

In order to search for a second 7r"^7r" kaon decay (-ft'2); all relevant tracks in the 
chamber — after removal of those originating from the decay already identified — 
are extrapolated to their points of closest approach to the K2 path computed 
from kinematics and xp. For each track candidate we compute the distance of 
closest approach to the K2 path. For each charge we take the tracks with smallest 
value of d as the K2 decay pions. We then determine the two track vertex and 
require that | m^+„- — rriK \ < 5 MeV and \ p^+^- — pk\ < 10 MeV/c. We also 

require — 50< -E^iss "Pmiss < 10 MeV^ and \J -E^iss + Pmiss < 10 MeV, where pmiss 
and -Emiss are the missing momentum and energy computed assuming A'2— ^tt+tt". 

More accurate values for the Ki^2 vertex positions, fi,2 and the collision point, 
fc, are obtained from a kinematical fit. The fit makes use of the constraint from 
the Kl and K2 directions ni 2 defined by: 

^ ^ (P+ + P-)l,2 - (P+ + P-)2,l + 

^'^ l(P+ + P-)l,2- (P+ + P_)2,l + P^r 

We then have 

fi = fc + /jiii, z = 1, 2 
and solve for fc and Zi 2 by maximizing the likelihood 

lnL= ^ lnPi(ri,fc + /ini) +lnPc(rc,fc), 

i=l,2 

where Pi and Pq are the probability density functions for rj and rc, respectively, 
as obtained from MC. The value of li^2 has been kept positive in the maximization. 
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The minimum of — In L is shown in Fig. 1 for data and MC. In order to maximize 




-InL 

Fig. 1. Distribution of the minimum of — In L of the kinematic fit for data and Monte 
Carlo for different decay channels. 

signal efficiency, improve At resolution, and minimize incoherent regeneration 
background, as a final selection requirement we retain events with — InL < 7.5. 

The time difference between Ki 7i~^tt~ and K2 tt+tt" is determined as 
At = \ti — t2\, where the proper time is ti = /j/(/3i7i) with /3j7j = p^Jrax- 

3.2 Determination of background 

The selected KsKl 7i^7i^7i^n~ events have a contamination of ~3.2% for 
At < 35rs dominated by regeneration on the beam pipe. Semileptonic Kl de- 
cays amount to 0.2% as determined from MC. Direct four pion production, 
e'^e~ — *• 7r"''7r~7r~''7r~, gives a ~0.3% contamination at the IP with At ~ 0, the 
region most sensitive to coherence loss. This contribution is obtained from the 

sidebands: 10 MeV < \/^^ii^7+"pLss < 20 MeV and - pk\ > I MeV. Ver- 

tex positions, total energy E^^,, and total momentum ^4,^ are used to distinguish 
between e'^e^^Tr'^ir^TT^iT^ events and semileptonic decays near the IP. We 
find 27±8 events, in agreement with the estimate of ~32 events from the cross 
section given in Ref. [21]. 

Incoherent (see Sec. 3.4) and coherent regeneration on the beam-pipe are included 
in the fit of the At distribution. The coherent regeneration amplitude, pcoh = 
|Pcoh|e*'^'=°'", is obtained from the time distribution for Kl{—^ Ks) — n^Ti~ decays 
after single K^s cross the beam pipe: 
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I{t)= (7r+7r-|/u(t) + Pcoh|A5(t)) 
= {n+n-\Ks)W\v-,-\'e~^-'- 



(8) 



Pcohl'e-^^* 



decay, using the same 
branching ratio [22]. 



PKr 



+ m\ 



K^s are identified by the reconstruction of a Ks- 
algorithm used for the measurement of the - 
In addition we require \pkl ~ Pn+n- I < 5 MeV and I-Ett+t 
mTT+TT- — '^T'kI < 5 MeV, where the Kl momentum, px^, is obtained from the 
Ks direction and p,^ and -Ett+tt- is the energy of the tt+tt" pair. These cuts are 
effective for the identification of KL—>-7r~^7r~ decays and coherent and incoherent 
/('2.~^-^5^7r"^7r~regeneration processes with a neghgible amount of background. 
Incoherent regeneration is rejected by requiring that the angle between p/^^ and 
p^+^- be smaller than 0.04. The residual contamination is ~3%, from the side- 
bands of the distribution in the above angle. Fitting the proper-decay-time dis- 
tribution, Eq. (8), we obtain |pcoh| = (6.5 ±2.2) x 10"^ and 0coh = (-1.05 ±0.25) 
rad. This result is stable against variations of the scattering angle cut and agrees 
with predictions [23] . Coherent regeneration in the inner pipe is negligible. Back- 
ground due to production of C-even neutral kaon pairs in two photon processes 
or (/o, flo) decays is also negligible [1,24-26]. 



3.3 Determination of the detection efficiency 



The overall detection efficiency is about 30%, and has contributions from the event 
reconstruction and event selection efficiencies. These efficiencies have been evalu- 
ated from MC. For the reconstruction efficiency, a correction obtained from data is 
applied. This correction is determined using an independent sample of KsKl 
TT+TT", TT/iz/ decays, tt/xz/ decays are identified by requiring \J -E^iss + Pmiss > 10 
MeV, Ip^I < 246 MeV and p\ + p*_ < 367 MeV, where p* is the momentum of 
the decay secondary in the kaon rest-frame, calculated assuming the vr mass hy- 
pothesis. We then compute the squared lepton mass mj-{mf+) in the hypothesis 
K 7r+/i-z/ {K vr-/i+z/), and require: 150 MeV^ < mf. + mf+ < 270 MeV^. 
The distribution of the time difference At = |ti — ^2! between the two kaon decays 
obtained for the KsK^ — > 7r"'"7r~, vr/xz/ sample is shown in Fig. 2, both for data 
and MC. The correction to the reconstruction efficiency from the MC is obtained 
from the data-MC ratio of the distributions in Fig. 2 and applied bin by bin as a 
function of At. 

In order to take into account the At resolution when fitting, a smearing matrix 
has been constructed from MC by filling a two-dimensional histogram with the 
"true" and reconstructed values of At. The efficiency correction and the smearing 
matrix are then used in the fit procedure as explained in the following section. 
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Fig. 2. At distribution for the KsKl — > 7r"'"7r~, vr/ii^ control sample for data (black 
points) and Monte Carlo (solid histogram). The expected background contamination 
from KsKl tt~^tt~ ,7rei' is also shown. The hatched area represents the Monte Carlo 
statistical uncertainty. 

3.4 Fit 

We fit tlie observed At distribution between and 35ts in intervals of At = ts- 
The fitting function is obtained from the 1(^1,^2) distribution given in Eq. (2) 
including the QM violating parameter (, or the QM and CPT violating parame- 
ters 7 and uj as discussed in Sec.l. To take coherent regeneration into account in 
Eq. (2), the time evolution of the single kaon is modified as follows: 

\Ks,L{ti)) = \Ks,L{ti)) + Pcoh\KL,s{ti)) (9) 

where pcoh is evaluated as explained above. We then integrate 1(^1,^2) over the 
sum ti + 12 for fixed At = |ti — ^2!; and over the bin- width of the data histogram: 

jAt 00 
/,(q)= J d{At) J I{ti,t2;q)d{ti+t2), 
(j-l)At At 

where q is the vector of the QM- and CPT-violating parameters. Finally, the 
observed At distribution is fitted with the following function: 

= iV I ^ s,, /,(q) | + N^'^I^^ + N,J,^,i (10) 



• Data 
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where rij is the expected number of events in the i bin of the histogram, Sij 
is the smearing matrix, and ej is the efficiency. A^, the number of KsK^ 
Ti^Tx~Ti^Ti~ events, and iV^*^^, the number of events due to incoherent regeneration, 
are free parameters in the fit. The time distribution fj^^ for the contribution 
from incoherent regeneration is evaluated from MC. The contribution from non- 
resonant e^e~ 71^71^71^71^ events is treated in a similar manner, except for 
that N^T^ is fixed to the value determined as in Sec. 3.2, rather than left free in 
the fit. The fit is performed by minimizing the least squares function: 



where iV^ is the number of events observed in the i bin and 6ei is the error 
on the efficiency, including the correction. Using Eq. (10) with the QM- and 
CPT-violating parameters fixed to zero. Am can be left as free parameter and 
evaluated. In this case, the fit gives 

Am = (5.61 ±0.33) x 10%"\ 

which is compatible with the more precise value given by the PDG [27]: 

Am = (5.290 ± 0.015) x 10^ s~\ 

For the determination of the QM- and CPT-violating parameters. Am is fixed to 
the PDG value in all subsequent fits. As an example, the fit of the At distribution 
used to determine (sl is shown in Fig. 3: the peak in the vicinity of At ~ 17 ts 
is due to coherent and incoherent regeneration on the spherical beam pipe. 

3.5 Systematic uncertainties 

As possible contributions to the systematic uncertainties on the QM- and CPT- 
violating parameters determined, we have considered the effects of data-MC dis- 
crepancies (in particular on the At resolution), dependences on cut values, and 
imperfect knowledge of backgrounds and other input parameters. The contribu- 
tions from each source to the systematic uncertainty on each parameter deter- 
mined are summarized in Tab. 1, and discussed in further detail in the following. 
Since the QM- and CPT- violating parameters are most sensitive to small val- 
ues of At, particular attention has been devoted to the evaluation of systematic 
effects in that region. The dependence of the detection efficiency on At is mostly 
due to the cut on — In T in the kinematic fit used to evaluate the vertex positions. 
We have varied the cut from 6.5 to 8.5, corresponding to a fractional variation 
in the efficiency of ~ 5%. The corresponding changes in the final results are con- 
sistent with statistical fiuctuations. For each physical parameter determined, we 
take the systematic error from this source to be half of the difference between 
the highest and lowest parameter values obtained as a result of this study. These 
contributions are listed in the first column of Tab. 1. 
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Fig. 3. At distribution from the fit used to determine (sL- The black points with 
errors are data and the soUd histogram is the fit result. The uncertainty arising from 
the efficiency correction is shown as the hatched area. 

Table 1 

Summary of systematic uncertainties 



Cut 


Resolution 


Inputs 


Coherent 


vr+vr vr+vr 


stability 






reg. 


bckgnd 


SCsL 0.007 


0.002 


0.001 


0.001 


0.020 


(5CooXlO^ 0.03 


0.01 


0.01 




0.11 


(57 X 10^1 GeV 0.4 


0.2 


0.1 


0.3 


1.3 


6 ^Lox 10^ 0.8 


0.1 


0.1 


0.3 


1.4 


(5$ja;xlO^ 0.4 


0.4 


0.2 


0.3 


0.1 



The QM- and CPT- violating parameters depend also on the At resolution. We 
have checked the reliability of the MC simulation on a sample of events with 
Ki=Ks K2=Kl by requiring I2 > 12 cm and removing the li > cut. Fig. 4 shows 
the Ks proper-time distribution for data and MC. From the negative tail of the Kg 
proper-time distribution we obtain the experimental resolution. We fit the data 
and MC distributions to an exponential function convoluted with the resolution. 
We obtain an rms spread of (1.152±0.020) ts for data and (1.1807±0.0036) ts for 
MC, i.e., the data and MC resolutions agree to within 1.4 standard deviations. In 
addition, we obtain a Ks lifetime, ts = (0.9030±0.056) 10"^'' s, in agreement with 
the world average value (0.8958 ± 0.0005) 10"^° s [27]. To estimate the resulting 
contributions to the systematic uncertainties on the QM- and CPT- violating 
parameters, the MC resolution is varied by ±5%, about three times the statistical 
uncertainty of the check. For each parameter determined, we take the systematic 
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error from this source to be half of the difference between the highest and lowest 
values obtained. These contributions are listed in the second column of Tab. 1. 
The third column gives the contributions of uncertainties on the known values of 
Am, r^,, and ?7+-, which have been propagated numerically. 

Contributions to the systematic uncertainties due to limited knowledge on |pcoh| 
and <^coh have been evaluated by varying the parameter values within their errors 
and are listed in the fourth column of Tab. 1. 

Finally the last column gives the contributions arising from the uncertainty on the 
level of background contamination from non-resonant e^e~ — ^ 7r+7r~7r"'"7r~ events, 
evaluated by varying the background parameter in the fits {N^.,^ in Eq. (10)) within 
its error. Note, however, that these contributions are included in the statistical 
uncertainties, rather than in the systematic uncertainties, in the statement of the 
final results. 

4 Results and conclusions 

From the fit we obtain the decoherence parameter values: 

CsL = 0.018 ± 0.040stat ± O-OOTgyst xV^of = 29.7/32; 

Coo = (O.IO ± 0.21stat ± 0.04syst) x 10"^ x^/doi = 29.6/32. 

which are consistent with C = cind no QM modification. Using the Neyman 
procedure [28], we derive the upper limits Csl < 0.098 and Coo < 0.50 x 10~^ at 




t/x 

s 

Fig. 4. Distribution of Kg proper-time distribution for data (black points) and Monte 
Carlo (solid histogram). 
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95% CL. Since decoherence in the K^K^ basis would result in the CP allowed 
KsKs — > TT+TT^TT'^Tr" decays, Eq. 4, the value for (^qo is naturally much smaller. 
In the model of Ref. [5], we find: 

A = (0.13±0.30stat±0.05syst)xl0-^^ GeV, A < 0.73x10"^^ GeV at 95% C.L.. 

All the above results are a considerable improvement on those obtained from 
CPLEAR data [4,5]. We have measured the 7 parameter: 

7 = (l.3l|^±0.4) X 10-21 GeV 

with xVdof = 33/32. From the above we find 7 < 6.4 x IQ-^i GeV at 95 % C.L.. 
This result is competitive with that obtained by CPLEAR [15] using single kaon 
beams. 

0.2 
0.15 

0.1 
0.05 



-0.05 
-0.1 

-0.1 -0.05 0.05 0.1 0.15 0.2 

Fig. 5. Contour plot of ^u; versus 3ft w at 68% and 95% C.L. 

The complex parameter uj has been measured for the first time. The result is 

= {l.ltli ± 0.9) X 10"^ Quj= {SAtt^o ± 0-6) X 10"^ 

with x^/dof = 29/31. The correlation coefficient between and ^Scj is 90%. 
Fig. 5 gives the 68% and 95% C.L. contours in the Qu, plane. The upper 
limit is |cj| < 2.1 X IQ-^ at 95% C.L. 

We do not find any evidence for QM or CPT violation. 
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